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ABSTRACT 

We calculate radio-to-X-ray light curves for afterglows caused by non-thermal emission from a highly 
relativistic blast wave, which is inferred from the gamma-ray flux detected in GRB 980425 and from 
the very bright radio emission detected in SN 1998bw. We find that the observed gamma-ray and radio 
light curves are roughly reproduced by the synchrotron emission from a relativistic fireball. The optical 
flux predicted for the non-thermal emission is well below that of the thermal emission observed for SN 
I998bw so that it will not be seen at least for a few years. The model predicts the X-ray flux just 
above the detection limit of it BeppoSAX for the epoch when it was pointed to the field of GRB980425. 
Therefore, the nondetection of X-ray and optical afterglows is consistent with the model. The models 
presented here are consistent with the physical association between SN I998bw and GRB980425, and 
lend further support to the idea that this object might correspond to an event similar to the "hypernova" 
or "collapsar" - events in which the collapse of a massive star forms a rotating black hole surrounded 
by a disk of the remnant stellar mantle. 

Subject headings: gamma-rays:bursts - supernovae:general - supernove: individual (SN I998bw) 



the star results in the formation of a system composed of 
a Kerr black hole and a disc of the remnant stellar mantle. 
Then, a required amount of energy may be extracted from 
the system by either neutrino annihilation or megneto- 
hydrodynamic effect, i.e., the Blandford-Znajek mecha- 
nism(Blandford & Znajek 1977), eventually producing a 
relativistic jet and a bulk mass ejection. In these scenar- 
ios, the high energy photons from the jet may be observed 
as a 7-ray burst and the bulk mass with a slower expan- 
sion speed would be seen as a supernova-like object like SN 
1998bw. Therefore, SN 1998bw and GRB980425 provide a 
good opportunity to test this class of models for gamma- 
ray bursters. In this Letter, we combine the observational 
data in different wavelength bands and try to see if they 
can be consistent with the hypernova or collapsar scenar- 
ios. We briefly give a description of a fireball model that 
we apply in §2. Then we compare the model prediction 
with the data in §3 and discuss possible interpretations 
and their implications in §4. Concluding remarks will be 
given in §5. 

2. BLAST WAVE MODEL 

We apply a simple fireball model for the blast wave dy- 
namics, as in Meszaros & Rees(1997) and Wijers, Rees, 
& Meszaros(1997). The blast wave with an initially large 
Lorentz factor T decelerates as it sweeps up external mat- 
ter, giving rise to non-thermal emissions at a decreasing 
characteristic frequency. The synchrotron radiation from 
electrons accelerated near the shock front is one of the 
most viable emission mechanisms, unless the density is 
too low and the energy transfer from protons to electrons 
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1. INTRODUCTION 

The discovery of the peculiar supernova SN 1998bw in 
the error box of the 7-ray burst GRB 980425 raised the 
possibility that at least some classes of 7-ray bursts may 
originate from supernovae(Galama et al. 1998, Kulkarni 
et al. 1998). The optical properties of SN 1998bw suggest 
that it was an explosion of a massive star with a kinetic 
energy 10 times larger than usual supernovae (Iwamoto et 
al. 1998, Woosley, Eastman, & Schmidt 1998). The radio 
light curves of SN 1998bw indicate the existence of a rela- 
tivistic blast wave associated with the supernova (Kulkarni 
et al. 1998). The 7-ray burst GRB 980425 was an unre- 
markable event in terms of spectral and temporal proper- 
ties. The time of occurrence of the burst coincides with 
that of supernova to within (+0.7, -2.0) days, and the 
supernova was found in the 8' error box of GRB 980425 
determined by BeppoSAX. No further evidence has been 
established so far to support the physical association of 
these two events, however, the chance probability of find- 
ing a supernova in the error box is estimated to be rela- 
tively small, ~ 9 x 10~ 5 (Galama et al. 1998). Provided 
that GRB 980425 was really associated with SN 1998bw, 
the energy radiated in 7-rays turns out to be ~ 8 x 10 47 
erg, assuming isotropic emission, which is 4 orders of mag- 
nitude smaller than in other bursts whose distances have 
been estimated so far(Galama et al. 1998). 

The possible supernova connection of the 7-ray burst re- 
minds us of "hypernova" (Paczyhski 1998) or "collapsar" 
(Woosley et al.1998; MacFadyen & Woosley 1998) scenar- 
ios for 7-ray bursters, in which the collapse of a very mas- 
sive star is involved. The scenario is that the collapse of 
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does not occur sufficiently quickly. The purpose of this 
Letter is not to give a unique and comprehensive model 
but to check the consistency of the physical association of 
GRB980425 and SN 1998bw. Therefore, we try a simple 
fireball model similar to those applied to GRB970228 in 
Wijers et al.(1997). 

Recently, Rees & Meszaros(1998) suggested the possi- 
bility that the slower parts of ejecta would catch up with 
the decelerated blast wave and may reenergize it. The en- 
ergy could be supplied to the blast wave during its prop- 
agation, resulting in afterglows that could be even more 
powerful than the 7-ray burst itself. In our analysis, we 
also assume a power-law T evolution so that T oc r~ n as 
in Rees & Meszaros(1998), which is realized in self-similar 
solutions (Blandford & McKce 1976). Since the blast wave 
radius r is related to the observer time t by r <~ 2ctT 2 , the 
radius and the Lorentz factor are given by 

rKt l/(2n+l) ) roct -n/(2„+l)^ ^ 

The comoving frame synchrotron intensity is given by 



I' 



/ a for v < v' and I' oc for v > v' , where v' 



is the break frequency v' m oc B'% m and a = 1/3, f3 = 
— (p — l)/2, and p is the index of the energy spectrum of 
nonthermal electrons N e ("f c ) oc "f^ p . The typical Lorentz 
factor of electrons is given by 7 c , m ~ {m p /m e ) ■ T oc T 
so that i/j oc T 3 . The equipartition between magnetic 
field and internal energies, B' 2 /8tt ~ r 2 n cx m p c 2 , yields 
B' « T. (Primes indicate values evaluated in the co- 
moving frame.) The observed flux at the break frequency 
is given by F Vm oc (ctL) 2 L 3 /^ oc i 2 r 5 /^ . The comov- 
ing intensity at the break frequency is given by I' Um oc 
<£ ,2 7c 2 mAr/(B'7 2 m ) oc n' c B'Ar oc L • L • rP" 1 = Tr? 

Therefore, the observed break frequency and the flux at 
the break frequency evolve as 



and 



rv' oc r 4 oc i-4«/(2«+D 



F v „ oci 2 r 6 roci-( 2 "- 3 )/( 2 " +1 \ 



(2) 



(3) 



respectively. 

The index n is determined by equating the ram pres- 
sures on the shocked region from the forward shock pf and 
from the reverse shock p ICV (Rees & Meszaros 1998). We 
assume, for simplicity, that the mass contained in the lay- 
ers outside the shell with T — Tf, M(T > Tf), is given 
by 

M(T > L f ) cx ry fl , (4) 
and a power-law density structure of external medium 

Pcxt oc r~*. (5) 

The mass that has caught up with the contact discon- 
tinuity when the blast wave reaches radius r is given by 



M r oc r ns . Then we have pf oc dM T /dr 



and 



p rov Poxtr 2 r 2 . In the latter, we assume the blast wave 
is adiabatic. From pf ~ p rcv , we find the index to be 



n = 



3-i 
s + 2' 



(6) 



For a uniform distribution of external matter ( t — 
) and a shell with a unique T ( s = 0), we have n = 
3/2, which reproduces the result obtained in Wijers et 
al.(1997): r oc t 1 / 4 ^ oc ^ 3/8 ,i^ m oc i" 3 / 2 ,^ oc t° = 
const. 

The monochromatic light curve is then given, if no 
beaming occurs, by 



Fu = F„ 



and 



oc t 



. -2n + 3 + 4na 
S = 2n + l ■ '<*" 



v V s > , _ -2n + 3 + 4n/3 
■ - - ; //,,, \ — j oc t , — -- 



F — F 



2n + 1 



, t ^> ti/) 

(7) 

where t v is the time of the light curve break at frequency 
v, which satisfies 



U> (u\ 1/2+1/4™ 



as seen from equations (2) and (7). 
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3. COMPARISON WITH THE DATA 

Table 1 lists the reports on the 7- and X-ray de- 
tections in the Wide-Field Camera(WFC) error box of 
GRB980425 (Pian et al. 1998a, 1998b). We converted 
these data to the fluxes in mJy (= 10~ 26 erg s _1 cm~ 2 
Hz -1 ) and plotted them in Figure 1 with open and filled 
circles, respectively. The observed optical light curve of SN 
1998bw in visual band(A - 5500 A) (Galama et al. 1998) 
and the radio radio light curve at 6 cm ( Kulkarni et al. 
1998) are also plotted with open and filled squares, respec- 
tively. A theoretical light curve of model C0138 (Iwamoto 
et al. 1998) for SN 1998bw, which is the thermal emis- 
sion from a subrelativistic ejecta heated by radioactivity, 
is shown with a solid line. Figure 1 also shows a set of 
afterglow light curves in X-rays(/i^ = 5keV) {dashed line), 
optical(visual band, which is at approximately A ~ 5,500 
A; dash-dotted line), and radio (A= 6 cm; dotted line). It is 
seen that the radio light curve agrees with the model pre- 
diction, although it shows a complicated behaviour before 
and around the peak. 

The gap of the peak fluxes between 7-ray and radio 
bands requires that F Vm should increase as a function of 
time. The slope of F Vm is chosen to be 1 /3 as shown in Fig- 
ure 1, which requires n = 1, thus s+t = 1. In order for the 
X-ray flux to fall below the detection limit at t = 10 4 93 s 
= 1 d, the X-ray light curve should decay faster than t -12 , 
which means [3 < — 1.15 if n = 1 as seen from equation (7). 
We choose (3 = —1.15 to draw the afterglow light curves 
in Figure 1 . Owing to the unknown evolution of random 
magnetic field strength B' or due to the other emission 
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processes involved such as self-absorptions and the inverse 
Compton scattering, the spectrum shape is uncertain and 
several breaks are even expected to appear in the spec- 
trum(Piran 1997). Therefore, we simply assume a = 
as adopted in Wijers et al.(1997) for our analysis in this 
Letter. 

From equation (8) with n = 1, one expects the time 
of the breaks at the X-ray, optical, and radio bands as 
t x = 10 - 9 %,i op ticai = 10 3 ' 5 i 7 ,t radio = 10 7 ' 2 t 7 , respec- 
tively. The flux of the optical afterglow is expected to be 
well below the observed optical flux of SN 1998bw; thus, 
it is no wonder that it has not been detected. The model 
predicts that the optical afterglow will not be seen at least 
for a few years, and then it will be below detection limits 
even with a largest scale of telescopes. In contrast, the ra- 
dio flux is still at a level close to its peak, and monitoring 
the manner of its further decay is crucial to distinguish 
the nature of the blast wave. The X-ray flux detected by 
BeppoSAX at the time of the burst is significantly large 
compared with the 7-ray flux itself so that it seems hard 
to explain. 

From the radio properties of SN 1998bw, Kulkarni et al. 
(1998) inferred that there exists a relativistic shock with 
r^l — 2 at t = 10 6 ' 5 s that is responsible for the radio 
emission. Equation (1) predicts T oc i -1 / 3 if we take n = 1; 
thus, we can estimate the initial bulk Lorentz factor Tq. 
The result is r ~ (1 - 2) x (10 6 5 /10 - 5 ) 1 / 3 = 100 - 200, 
where we set the time of the burst as 10 0,5 s. Kulka- 
rni et al. (1998) also gave an estimate of the mass of 
the relativistic shock M e j ~ 1O~ 5 M . Then the initial 
energy of the blast wave is estimated to be less than 
r M cj c 2 - (1.8 - 3.6) x 10 51 erg. 

4. DISCUSSION 

Neither of the X-ray transient sources detected by 
Narrow-Field Instruments (NFI) on BeppoSAX coincides 
with SN 1998bw in their positions. One of them (1SAX 
J1935. 0-5248) was reported to have a constant flux, but 
the other one (1SAX J1935.3-5252) has been fading(Pian 
et al. 1998a, 1998b). If the latter was indeed an afterglow 
of GRB 980425 and the burst occurred at the position, 
the fireball model predicts an optical afterglow and a sub- 
sequent radio afterglow with a detectable level of fluxes. 
Since the positions of the two X-ray sources were 3' away 
from that of SN 1998bw, they would have been visible 
there. The lack of such optical and radio afterglows may 
also indicate that the sources had nothing to do with GRB 
980425. 

Kulkarni et al.(1998) reported the radio light curves in 
four different wavelength bands (Figure 2). If we closely 
look at the light curves, we see that they show compli- 
cated structures before and around their maximum epoch, 
but rather clean spectrum and temporal evolutions later. 
The irregularities seen in the early light curves may pos- 
sibly be the result of inhomogencity in external matter, 
absorption by dust in the host galaxy, and radiative pro- 
cesses other than synchrotron radiation. The light curves 



are well fitted by power-law decays with the same expo- 
nent, 6' = —1.67. On the other hand, the flux ratios 
between different wavelength bands give a spectral index 
(3 = —0.7, which is close to the values reported by Kulkarni 
et al.(1998). These requires n = 3.2, which does not agree 
with the value we chose to fit the light curves in Figure 1. 
Such a large value of n leads to a negative t, which may in- 
deed suggest a density inhomogeneity in external matter, 
although it is necessary to do more detailed modeling of 
the radio light curves with a realistic treatment of hydro- 
dynamics including other potentially important radiative 
processes (Nakamura et al. 1998). 

5. CONCLUSIONS 

We examined the optical and radio light curves of SN 
1998bw and 7- and X-ray fluxes observed in GRB980425. 
Under the hypothesis that the two events are of the same 
origin, we compare the light curves that a simple fireball 
model predicts with these observations. As a result, we 
find the following interesting facts that favor the possibil- 
ity that GRB 980425 was physically related to SN 1998bw. 

1. The radio emission from SN 1998bw and the 7-rays 
(BATSE, BeppoSAX) from GRB980425 can be interpreted 
as a single event based on a simple fireball model, although 
the X-ray flux reported by BeppoSAX team is a bit too 
large and hard to be reconciled with a model as simple as 
the one adopted here. 

2. The X-ray flux predicted by the model is only 
marginally detectable at the epochs BeppoSAX NFI was 
pointed to the field, which is consistent with the nonde- 
tection of any X-ray afterglow in GRB980425. The pre- 
dicted X-ray flux falls below the detection limit of ASCA 
in ~ 2 — 3 months; therefore, it is not likely that we can 
detect decaying X-ray emissions from the burst by further 
observations. 

3. Observations suggest that there seems to be a relativis- 
tic blast wave and a subrclativistic bulk mass ejection in SN 
1998bw. The former might have caused the 7-ray emission 
from GRB980425 and the radio emission from SN 1998bw, 
while the latter corresponds to the optical emission from 
SN 1998bw. This picture is consistent with the hypernova 
or collapsar scenarios for 7-ray bursters, which involve the 
collapse of a massive star. 

Although the arguments here are not strong enough to 
claim the physical association between SN 1998bw and 
GRB 980425, it is worth noticing that the relativistic shock 
required to explain the radio emission from SN 1998bw 
could be interpreted as a decelerated "fireball" that ini- 
tially had a higher bulk Lorentz factor and was able to 
radiate 7-rays observed in GRB980425. 
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useful comments to improve the manuscript and the fig- 
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Note.l: The late-time optical light curve of model C0138 is calculated by extending the exponential decay of the light 
curve of the same model as presented in Iwamoto et al. (1998). The previous version of the light curve was not accurate 
enough, and has been replaced by a correct one. 



Note. 2: The positions of the two NFI X-ray sources have been revised (GCN Circ. 155, which says that "one of the NFI 
sources was at 50" from SN 1998bw and therefore consistent with it"). A following GCN Circular (No. 158) reported that 
the flux from the 'consistent source' has shown a moderate decay ( approximately by a factor of two from April-May to 
November in 1998). The slow decline could be a possible sign of circumstellar interaction in SN 1998bw(Note that the 
last statement is merely our guess yet). 



Time 


Energy Band (keV) 


Flux (or Fluence) 


Source 


0-25 s 


24-1,820 


(5.5 ±0.7) x 10 4fci erg s" 1 


CGRO BATSE 1 


0-30 s 


2-28 


3 Crab 


BeppoSAX WFC No.2 2 


~ 1 d 


2-10 


- (1.6 ± 0.3) x 10~ 13 erg s" 1 cur 2 


BeppoSAX NFI 2 - 3 


- 1.92 d 


2-10 


< 1.2 x 10~ 13 erg s^ 1 cm~ 2 


BeppoSAX NFI 2 ' 3 


- 6 d 


2-10 


< 1.0 x 10~ 13 erg s _1 cm~ 2 


BeppoSAX NFI 2 ' 3 



Table 1 

Voted from Galama et al. (1998), 2 Pian et al. (1998), 3 1SAX J1935.3-5252 
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Fig. 1. — Light curves in 7-ray, X-ray, optical, and radio wavelength bands predicted by a fireball model. The detection 
limits of BeppoSAX and ASCA satellites are indicated with bold dashed lines. The optical light curve of SN 1998bw and 
its model C0138(Iwamoto et al. 1998, see Note.l) are also plotted. 
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Fig. 2. — Radio light curves of SN 1998bw and at wavelengths, 3cm(filled squares), 6cm(open squares), 13cm(fillcd 
circles), and 20 cm(open circles) (Kulkarni et al. 1998). The early parts of the light curves show irregularities in 
spectral and temporal evolutions. However, the later decaying parts are well reproduced by power law decay curves with 
5' = — 1.67,(3 = —0.7. Solid, dotted, short-dashed, and long-dashed lines correspond to 20, 13, 6, and 3 cm bands. 



